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We investigate a class of matrix model which describes the dynamics of identical par- 
ticles in even dimentional space. We show that the degrees of freedom after some 
constraints are implimented is proportional to particle number and consist of those for 
positions and internal degrees. The particle dynamics is given by the metric on the 
smooth moduli space. The moduli space metric for two particles is found. The size 
of tightly packed particles grows like \/iV. Our matrix model is related to the ma- 
trix model for fractional quantum Hall effect, the ADHM formalism of U{\) instantons 
on noncommutative space, and supersymmetric DO branes on D6 branes with nonzero 
B-field in type IIA theory. 
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Recently, there have been some studies to understand the dynamics of identical particles in terms 
of the matrix gauge theory. For N particles in D + 1 dimensional space, the particle coordinates 
become D + 1 square matrices of size A^. The particle exchange symmetry becomes a part of the 
gauge symmetry and so is included in the dynamics classically. Usually the matrix theory has more 
degrees of freedom than ones necessary for particle dynamics and further constraints are imposed 
to reduce the degrees of freedom. A typical example is the dynamics of instanton in terms of 
the matrices in the ADHM construction [Q]. These matrices are constrained by the three D-term 
constraints and by the Gauss law and so the number of net degrees of freedom is proportional to 
the instanton number |^, |2[. The number of degrees of freedom for a single instanton is the sum 
of that for the position and that for the internal degrees. Another example is the Chern-Simons 
inspired matrix model which describes the fractional quantum Hall fluid. ^ . 

In this work we study a matrix theory which describes identical particles in even dimensional 
space. Our matrix theory is quadratic in time derivative and so naturally describes the nonrela- 
tivistic Newtonian dynamics of identical particles. We count the number of degrees of freedom left 
after the constraints are imposed, which turns out to be proportional to particle number. For the 
relative moduli space for two particle case, we extend our result on four dimensions to other even 
dimensions. The size of tightly packed N particles grows like \/]V. 

Our theories appear in many context. On two dimensional space it is somewhat related to the 
first order matrix theory introduced by Susskind and Polychronakos. On four dimensional space 
with single flavor it is related to the ADHM formalism of U{1) instantons on noncommutative 
four space. On six dimensional space it describes the low energy dynamics of super symmetric DO 
branes on D6 branes, which is possible with nontrivial background NS B field in type HA theory 
as discussed by Witten . Some of issues studied here has been discussed in somewhat simplier 
level in Ref. g. 

The ADHM description of the low energy dynamics of instanton is a gauge theory of matrices 
with additional constraints, which can be solved in terms of the moduli parameter. Once the 
constraints are solved and the gauge degrees of freedom are removed, the instanton moduli space 
metric describes the motion of instantons. The metric on the moduli space is singular and the 
instanton dynamics is incomplete. 

The ADHM formalism arises naturally from DO brane on D4 brane point of view. Its low energy 
dynamics of DO brane is given by the Yang- Mills matrix mechanics model with eight supersymme- 
tries. The relevent fields are adjoint and fundamental scalar fields. These matrices are constrained 
in the ADHM formalism by the three D-term potentials 0. 

When nonzero background B field introduced on D4 branes, their field theory becomes Yang-Mills 
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theory on noncommutative space [Q, ^, which can be obtained from a certain hmit of open string 
theory on D4 branes with background B-field. In this case, one can have nonsingular moduh space 
for instantons Q. (See Ref. ||l^, 11 1 for a recent review.) The B-field appears as the FI term in the 
matrix theory, which blows up the singularities of the instanton moduli space with a finite scale. 



The ADHM formalism is simplified when the gauge group is U{1) |12| , 13]. Our class of the matrix 
model includes only the ADHM formalism for the U{1) group. 

In two dimensional space Susskind and Polychronakos proposed and studied the matrix model 
for charged particles on the lowest Landau level and so their Lagrangian is first order in time 
derivative and found to be related to the noncommutative Chern-Simons theory ||,||. The space 
of matrices become a classical phase space. After the Gauss law is imposed and the gauge degrees are 
removed, the remaining degrees of freedom are proportional to the number of particles. Especially 
the quantum mechanical eigenvalues and eigenfunctions of a quadratic Hamiltonian can be found 



exactly |14 



On the other hand our matrix model is second order in time derivative and so it describes the real 
time evolution of particles. After reducing the degrees of freedom, the motion of identical particles 
is described by the coordinates of the moduli space. The moduli space, which is the coordinates 
space in our model, appears as the phase space in the matrix model of Susskind and Polychronakos. 

In the type llA theory, DO branes feel repulsive force from the parallel D6 branes. However when 
there exists a uniform background B field above critical strength, DO branes can be attracted to D6 
branes and form a BPS configuration. The DO brane physics is described by the matrix model 
This is exactly the matrix model we study in the six dimensions. 

Noninteracting identical particles on fiat have a singular configuration space {R^)^ / Sn where 
Sn is the permutation group of N particles. The configuration space is singular when particles 
come together. On even dimensional space our matrix model is a natural blow-up procedure of 
these singularities. There is a natural length scale for this blow up, making particles to carry 
effectively finite size core. When many identical particles described by the matrix come together, 
their kinetic energy, which determines the moduli space metric, makes them to behave somewhat 
like incompressible fiuid even though there is no force between them. In our case the total volume 
occupied by tightly packed particles in D = 2k dimensions grows as A^'^, which is faster than 
the incompressible gas except k = 1. As we will see, the kinetic energy also entails the long range 
interaction between particles. Thus, one could regard the matrix model to describe a peculiar class 
of particles living on even dimensional space. 

We start with the matrix model description applicable for all these cases. The matrix model for 
A^ particles on 2k dimensional space is described by A^ x A^ complex matrices Zi with i = 1, ...,k 
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and N X M complex matrix ip. The kinetic term for these matrices is 



£ = tr {DoZiDoZ, + Dq^jDo^) . (1) 

The U{N) local gauge transformation leads to Zi — > UZiU, tp Utp and Aq UAqU — idoUU. 
The U{M) global flavor symmetry leads to ■0 ^ The matrices Zi,tp have 2kN'^ + 2NM 

real parameters. Out of which A^^ are gauge parameters, leading to {2k — 1)N'^ + 2NM physical 
parameters, which is clearly much larger than the number of particles. The spatial rotational group 
S0{2k) is broken to U{k). (See for a recent investigation of the model with k = 1/2 and so particles 
living on a line where the matrix degrees of freedom are proportional to the particle number |15].) 

Additional constraints are needed for these matrices. We choose them to be 

k 

J2[Zi,Zi] + i^i^ = CiN, (2) 

2=1 

[Zi,Zj]=0 i,j = l,2,...,fc., (3) 

These can be regarded as the minimum of a potential of matrices. In the theory with four super- 
symmetries the first one (^) can be regarded as the minimum condition of the D term potential 
and the second one as the minimum condition of the F term potential. 

For the two dimensional case, there is only one constraint (^). The number of the maximal 
supersymmetry is four and the constraint is the minimization condition of the D-term potential. 
For the four dimensional case with a single flavor M = 1, the above constrains are reduced ver- 
sion of the ADHM constraints for the U{1) instantons on noncommutative four space, and so the 
number of underlying supersymmetry is eight. For the theories of the six dimensional system, the 
number of maximal supersymmetry is again four and the constraint (^) comes from the superpo- 
tential ti Zi[Z2, Z^]. For other cases there is no obvious supersymmetry is as there is no suitable 
superpotential for the constraint 

For one particle with = 1 Eqs @ and @ have the trivial solution where Zi are arbitrary 
numbers and 

V'V^ = CliV. (4) 

The space is 5*^*^""^ and becomes CP^'^^^ after we mod out by gauge symmetry U{1). The 
particle position Zi has 2k parameter and the internal degrees of freedom has 2(M — 1) parameters. 
The total degrees of freedom for a single particle is then 2k + 2M — 2. 

The super symmetric quantum dynamics of particles with M > 2 is possible on six dimensions. 
The quantum mechanics of the sigma model with moduli space metric has 4 supersymmetry as the 
maximal supersymmetry. For a single particle with nontrivial internal degrees of freedom, there 
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are M quantum mechanical ground states because they are normahzable harmonic forms |16| and 



their total number can be identified with the Euler number of CP^^ ^ by Hodge theorem [p!7| ]. 

When we consider many particles N > 2, one may wonder how degrees of freedom are there. 
Naively we expect that each particle carries the same number and so the total number of free 
parameters is the sum of that for constituent particles. However it is not obvious at all from Eqs. (§) 
and ^ imposed on the matrices Zi,ip. Eq.(|3|) seems over constraining. Indeed the constraints in 
Eq. dll) are not independent. In Appendix A we show that the correct counting leads to the 
N{2k + 2M — 2) degrees of freedom as expected. The original kinetic term then defines a Kahler 
space of dimension N{2k + 2M — 2) with the induced Kahler form from the flat Kahler form 
corresponding to the kinetic energy. Our matrix model is a natural method to define a class of 
Kahler spaces, which is related to the blow-up of the configuration space of identical particles. 

It is very hard to solve the constraints for arbitrary number of particles and the internal 
degrees of freedom characterized by M — 1. When M > 2, there are some internal degrees of 
freedom and so one can overlap the particles at a same point. Especially with M > N, one can 
solve Eqs. (|2|) and (^) with Zj = and the ip space satisfying Eq. (Q) becomes Grassmannian 
manifold U(M)/U{N) x U{M — N) after moding out the gauge group, as identified by Witten [^. 
When there is no internal degrees of freedom, N particles act like having a 'hard core', given not by 
force but by the metric of the kinetic term. When there are enough internal space M > 2, particles 
of number less than or equal to M can come together to the top on each other. 

On two dimensional space k = 1 with M = 1, the most general solution of the constraint (^) 
has been known before. (See for example Ref. [^.) A further discussion of our model on two 
dimensions and its relation to Polychronakos system is given in Appendix B. In four dimensional 
space the solution for two particle case has been found before [18| and has been used by us to show 



the moduli space metric to be the Eguchi-Hanson metric. In the following we extend this analysis 
to arbitrary even dimensions. 

With no internal degrees of freedom M = 1 and so there is only single flavor ip. In this case 
one can find the generic two particle solutions {N = 1) for these constraints by triangularizing the 
commuting complex matrices Zj simultaneously by unitary transformations. 



= + v. = (5) 

where 

Q= ^»l^^l >0, b= ^ . (6) 

Notice that the remnant f7(l)^ gauge is used for fixing the phase of a and b. We can identify wi 
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and Zi as the gauge invariant position of two particles at all separation of two particles not only 
at large separation since Wi it Zi/2 are the eigenvalues of Zi. There is still the remaining discrete 
gauge symmetry 

/ —b v2ab \ 



'2ah b 

which changes the sign of Zj, leaving ^ invariant. This discrete gauge symmetry identifies the 
positions of two identical particles and making the relative moduli space at large separation to be 
i?^^/Z2. The orbifold singularity at the origin gets blow-up in the matrix theory. At the coincident 
limit Zi = 0, the space of the solution (|^ becomes a submanifold S'^^~^. We have to mode out this 
by the U(l) gauge transformation diag(e*^, 1) because this transformation gets restored at this limit 
\/l — b = 0. Thus, the submanifold of moduli space at the coincident limit is S'^^~^ /U{1) which 
can be identified with CP^~^ because S'^'^"^ is U{\) fiber bundle over CP^~^ as we know [^]. Here 
CP^ denotes just one point. 

To obtain the tangent vectors on the moduli space, we consider the infinitesimal variation of the 
above solution made of the infinitesimal change in moduli parameters and the infinitesimal gauge 
transformation, 

6Zi = dZi — i[5a, Zi], 6ijj = dip — i6aip . (8) 

We need to fix the gauge variation of the above solutions under small variation, which can be 
achieved by solving the background gauge condition. The initial configuration for the matrices is 
the solution (|5[) and its 'initial velocity', which is characterized by the variation. It should satisfy 
the Gauss law constraint, which is identical to the background gauge fixing condition, 

k _ 

Y,{[6Zi, Zi] - [Zi, 5Z^]) + 5i>i^ -i;5i^ = {). (9) 

i=l 

This condition fixes the infinitesimal gauge transformation uniquely, 

• 7 / da— da \/2da 



where da = Zidzi/ (2C) and da = Zidzi/{2(). Knowing the tangent vector {6Zi,6ip), we can obtain 
the metric in the usual way 



ds^ = tr 6Zi6Zi + dil^S^^ . 



(11) 

Incidentally, Sip = in our two particle case. Then, the moduli space metric becomes 

ds^ = dsl^ + dsl^i, (12) 
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where the center of mass and relative metrices are 



dScm — 2^ dWidWi, (13) 



i=l 



"^'rel = ^2 dZidz^ - ^ Z.Z.-dZidz,- (14) 



with = I-Sjp- At the large r = \/J2i \ hmit the metric approaches that of the flat space 
with the correction of order 1/r^, which is a sign of long-range interaction. The above metric ds^^i 
for any /c is a Kahler metric as gjj = didjIC with the Kahler potential, 



^.. = iv/;^+^lnf!±g!-^^. (15) 

The Kahler form could be obtained by -fC = ^ddJC. This metric has f/(A;) holomorphic isometry 
transforming Zi — > UijZj, which is a subgroup of the spatial rotation S0{2k). In fact, the metric 
of the k = 2 case has hyperKahler structure and is identified with the standard Eguchi-Hanson 
manifold |18], where the isometry SU{2) C U{2) becomes triholomorphic. Basically the A: = 1,3 



cases are related to super symmetric theory of four real supercharges but the k = 2 case is related 
to theory of eight real supercharges. Therefore, the metric of the k = 2 case has the additional 
structure and is Ricci flat. Note that the metrics of the A: = 1, 3 cases are not non Ricci flat Kahler. 

To explore the region near the origin of the relative metric r ~ where two particles come 
together, we first note that the Fubini-Study metric on CP^^^ can be given by 

dsps = -2 ( X! dzidzi - ^ X! ZiZjdzidzjj . (16) 
After the following transformation Zi = Xqi, i = 1,2, k — 1, and Zk = X with p = ( X]f=i kiP)^''^! 



the above metric (16) becomes the standard one on CP 



k-l 



,2 _Ei=idqidqi EijliqiQjdqidqj 

"^FS — 2 2\2 ■ ^ ' 

The metric ds'j^'g is normalized so that the Ricci tensor satisfies R^j = 2k5fj. The CP^~^ space is 
Kahler with Kahler form K, which is locally exact, K, = dA. On E?^ , J^i^Zidzi = dr"^ + r'^dU'^k-i 
with ^2k-i being the metric on the unit S'^'^"^. The unit S'^^~^ sphere can be identified with U{1) 



fiber bundle over CP and its metric becomes [1£] 



dnlk-i = {dO + Af + dsls- (18) 
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Thus, the relative metric (14) becomes 



where the range of 9 is [0, vr] instead of [0, 27r] as we have identified Zi with — Zj. With change of 
variable v = r"^ the above metric becomes ds^ (df^ + u^d(20)^)/(8C) + 2(ds'pg near r 0. Thus, 
the above metric becomes the smooth metric of x CP^~^ near the origin. 



By the coordinate transformation 



^' = J+c^ (20) 



we can put the metric and Kahler potential in the following form 



dsli = -^+u\dsh + {l-^,){d9 + Af), (21) 
1 — " 



ICrel = n2 + ^ln(^^). (22) 



For the k = 2 case there exists unique self-dual middle harmonic form which corresponds to 
the supersymmetric normalizable ground state of the maximal eight super symmetric extension |l^ 
provides the wave functions of the threshold bound state of two U{1) instantons on noncommu- 
tative four space. For the k = 1,3 cases, it is not clear whether there exists any normalizable 
harmonic forms which correspond the supersymmetric normalizable ground states of the maximal 
four supersymmetric extension. 

Having studied two particle case, let us consider briefly particle case with M = 1. As all Zi are 
commuting due to the constraint (^), we can put all Zi as upper triangular matrices simultaneously. 
From the analogue with two particle system, we expect that all diagonal elements of Zi vanishes 
when all particles are packed tightly together. (Here we put the center of the mass position at 
the origin.) Then we can find the solution of two constraints (|2|) and (^). Only nonvanishing 
components of Zi are 

{Zi)a,a+1 = V^CUai a = 1,2, N - 1, (23) 
V'TV = (24) 

where Ua is /c-dimensional unit complex vector. We have used one of the diagonal t/(l)'s subgroups 
to make "0 vector real. Further use of the diagonal ?7(1) subgroups leads to identification Ua ~ 
g«(ea-9a+i)^^^ making each Ua vectors to belong to CP^"^ . Thus the solution space of N particles 



7 



for the relative motion becomes {CP ^) ^ when all particles come together. For this solution 

- N(N-l) , , 

tr Z^Zi =Y.<= <■ (25) 

a=l ^ 

While the particles are indistingushable, we could regard the radius square of the a-th particle is 
(a — The radius of the outmost particle is ^(A^ — 1)C and so the total volume in 2fc dimension 
will be proportional to N^(^ in large N limit. For particles in two dimensions, the minimum 
volume grows like which shows that they are incompressible. For particles in higher dimensions 
{p = 2k), the volume grows faster than N, showing the existence of additional mechanism. It 
shows nontrivial interaction between particles in short distance. Such phenomena is not new: the 
volme of nonabelian core of tightly packed N BPS monopoles grows as in three dimensions. 
The solution space of the packed particles, {CP^~^)^~^ is trivial for the two dimensional case and 
so the configuration is rotationally invariant. That is nontrivial for the higher dimensional case, 
and so is not invariant under U (k) subgroup of the rotational group. 

In summary we have explored a class of matrix model which describes identical particles in even 
dimensions. The number of free parameters of the theory turns out to be the sum of that for 
individual particles. The dynamics is described by the moduli space metric on the blow-up space 
of the singular space (C^)^ /Sn- We studied the two particle moduli space in detail. We show that 
the minimum volume of A'' particles grows like N^. 

Our model in two spatial dimensions and in four dimension with larger number of flavor M > 2 
docs not have the direct interpretation as the solitons in the field theory or the string theory. It 
would be interesting to find such field theoretic model. Finally, it would be a challenge to find 
similar matrix model for supersymmetric DO branes on D8 branes with nonzero B field. 
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Appendix A. 



To find out the number of free parameters in Eqs. (|2|) and (|^), let us start by noting all commuting 
square matrices can be triangularized simultaneously by unitary matrices. We triangularize all Zi 
matrices. After triangularization, we count the number of free real parameters. To do that, we want 
to show that [Zj, Zj+i] = for i = 1, .., /c — 1, then [Zj, Zj\ = for arbitrary i, j. Let us first start 
by showing that this is true for three upper triangular complex matrices ^1,^2,^3. 

Lemma. [Zi, Z2] = Z-^ = implies [Z3, Zi] = for generic upper triangular matrices. 

Proof) We prove this by induction. Suppose the result holds for kxk matrices. Then for (A; + 1) x 
[k + 1) Zj's, we write 

ry _ ( Si Vi\ 

^'-\ z J ' 



(A.l) 



where Si is a k x k upper triangular matrix, Vi is a k dimensional column vector and Zi is a number. 
Then 

[Zi Zj] = ^ '^■'^ ^^^^ ^-^^^ ~ ^^^^ ~ ^^'"^ ^ (A 2) 

Now from [Zi, Z2] = [Z2, Z^\ = we have [Si, S2] = [5*2, ^3] = 0. Then by induction hypothesis, it 
implies [^3, 5i] = 0. So the only thing need to prove is 

S-iVi + Z1V3 - S1V3 - Z3V1 = 0. (A.3) 

Indeed, this can be easily shown using the corresponding equations from (1,2) and (2,3) pairs. More 
precisely, after short calculations we get 

{S2 - Z2)iS3Vi + Z1V3 - S1V3 - Z3V1) = 0. (A. 4) 

But in general the determinant of 5*2 — Z2 is not zero since det(52 — Z2) = Yla=iiiS2)aa — Z2) and 
the diagonal components of 5*2 are not constrained by the commutators. Q.E.D. 

For four triangular matrices such that [Zi, Z2] = [Z2, Z3] = [Z3, Z4] = 0, the above lemma shows 
that [ZijZ^] = 0. Thus by using the above lemma again for ^1,^3,^4, we see [^1,^4] = 0. By 
induction, then it should be true for arbitrary number of Zi's. 

For each triangular Zi, there are A^(A^ + 1) real parameters. The total number of parameters of 
triangular Z'^s and are kN{N + 1)+2NM . The remaining gauge parameter after triangularization 
is and Eq.(^) imposes N"^ constraints. Due to the above argument there are {k — 1)N(N — 1) 
constraints on [Zi, Zj+i] = 0. Then total free parameters are 

kN{N + 1) + 2AM - - N-{k- 1)N{N - 1) = N{2k + 2M - 2). (A.5) 
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Appendix B. 

Here we explore the two dimensional case in detail. For k = 1, M = 1, the general solutions of 
the constraint @ for any have been found by introducing two Hermitian matrices X, Y such 
that Z = {X + iY)/\/2. After diagonalizing one of them by the unitary transformation, we use the 
remnant U{1)^ gauge freedom to put the column vector ip to be real positive. The solution of the 
constraint (0) is given as 



Xij 



) Xi for i = j 

[ -K/ivi - Vj) for ^ / J ' 
Yij = yAj, il^i = \/C. (B.l) 

In this coordinate the rotational symmetry is not manifest. (While the diagonal components of the 
triangularized Z matrices is rotationally manifest and the distance between particles are obvious, it 
is much hard to solve the constraint (§) in this coordinate.) Since the general solution is obtained, 
we might say that the moduli space metric can be obtained in principle and the multiparticle 
dynamics could be understood. However it turns out that it is technically somewhat difficult. Here 
we consider only the two particle case and find its explicit moduli space metric in this form of the 
solution. Again the Gauss law determines the infinitesimal gauge parameter 5a uniquely. After 
some calculation, the metric for the relative motion is given by 



2 I 2 

X + y 



dx"^ + 



1.2, 4C\2l. 2 SC^X 



2 I ^ / 2 I \2 

y y 



V 

where x = xi — xi-, y = yi — y2- (In this case it turns out Sip = 0, too.) 



dy"^ - ^^^dxdy ) , (B.2) 



Although {xi,yi) can be identified as the position of particles at large separation, its meaning at 
short distance changes. We can relate {x, y) to the relative coordinate z by taking the eigenvalues 
oi Z = {X + lY) I One can see that these two coordinates (x, y) for moduli space metric do not 
have one to one correspondence though one form of the metric can be obtained from the other by 
the direct coordinate relations. 



We need to elucidate the Polychronakos system at this point because Polychronakos has studied 
the first order time derivative Lagrangian Q 

£i = Ktr (iZDoZ + i^jDotp - C^o ) - w^tr ZZ, (B.3) 

where DqZ = Z — i[AQ,Z], and DqiP = tp — iA^ip.. The auxiliary field leads to the Gauss 
law constraint [Z,Z] + tp'ip = C^n, which is identical to the constraint Classically there are 
2N'^ + 2N degrees of freedom in Z and tp. Since there is only a first order term in time derivative, 
the space of Z, ip forms a phase space. A^^ parameters are constrained by the above Gauss law 
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and another N"^ parameters can be absorbed to the gauge parameters, leaving free 2N parameters 
in the phase space. These 2N parameters denote the N particle location on two dimensional 
noncommutative plane. The eigenvalues and eigenstates of the Hamiltonian H = up'iv (ZZ) are 



completely solved |14]. 



The first order kinetic energy (|B.3| ) leads to one form A = ti {ZdZ + iljdip) on the 2N dimensional 
moduli space, making it to be a phase space. As dA = tr (SZ A 5Z + 5^ A dip) regardless of 5a, the 
Kahler two form of our moduli space gives the symplectic structure of the Polychronakos system. 
Our moduli space becomes the phase space of the first order system. It would be interesting to 
work out the quantum problem in terms of the moduli space coordinate. 
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